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Abstract 

Elastic properties of Ti based /3-alloy were studied by the method of the model structure first prin- 
ciple calculations. Concentrational dependence of Young modulus for the binary /3-alloy Tii-^Nb^, was 
discovered. It is shown that peculiarities visible at 0.15-0.18% a.u. concentrations can be related to the 
different Nb atoms distribution. Detailed comparison of the calculation results with the measurement 
results was done. Young modulus for the set of the ordered structures with different Nb atoms loca- 
tion, which simulate triple /3-alloys Tio.518 Zro.297Nbo.i85 and Tio.297Zro.5i8Nbo.i85 have been calculated. 
The results of these calculations allowed us to suggest the concentration region for single-phase ternary 
/3-phase alloys possessing low values of Young's modulus. 

Introduction 

Today the /3-titanium alloys are widely used for medical applications due to 
their high bio-compatibility and untoxicity. The elastic moduli of these alloys 
have to be rather low. Elastic properties of the /3-alloys of Tii_^Nb^ system have 
been studied extensively. It is known that at Nb content below x=0.27 these 
alloys are not in a single-phase condition and contain martensite phases p-J. At 
higher Nb content single-phase conditions can be fixed. Martensite phases usually 
lead to a considerable increase in elastic constants of these alloys. Elastic moduli 
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of the single /3-phase apparently decrease with lowering Nb content, as the /3- 
phase turns less stable. Alloys with low elastic moduli of about 30-40 GPa are 
promising materials for implants. So, single-phase /3-alloys with Nb content below 
0.27 could be such material, however, as mentioned above they are practically 
impossible to obtain, and there is no experimental data on the /3-phase elastic 
moduli dependence on Nb content in Tii_^Nb^ alloys at x<0.27. It is worth to note 
that knowing this dependence is an important issue since it concerns changes in the 
elastic moduli of binary Tii^^iNb^; alloy when Ti atoms are partially substituted 
by other atoms. In this regard single-phase Tii_3;_j^ZryNbj; ternary alloys with 
low elastic moduli are of interest. In contrast to Tii_2,.Nba; alloys, the ternary 
alloys can be fixed in a single-phase condition at considerably lower Nb contents 
(up to x=0.14) [2j. It is well known also that Tii_^._yZryNb^. alloys have lower 
elastic moduli as compared to Tii_^Nb^ alloys. Therefore, partial substitution of 
Ti by Zr in Tii_^Nb^ alloy is expected to have a double effect - it will lower Nb 
content necessary to fix a single-phase state, and it will reduce the elastic moduli 
of the ternary alloy. At the first stage of search for Tii_j._yZryNbj; alloy optimal 
compositions, the information on the elastic moduli of binary Tii_^Nb^ /3-alloy in 
a wide range of Nb content should be obtained. Such information is valuable also 
for low Nb contents at which the alloy can not be practically fixed in a single- 
phase condition. In this regard first-principles calculations of electronic structure 
and elastic properties of the alloy are useful. High-precision quantum-mechanical 
calculations allow to study in detail various factors affecting the elastic moduli of a 
material. Apparently the most important factors are Nb content and distribution 
of Nb atoms in the BCC lattice. Whether the elastic moduli dependence on Nb 
content represents a monotonic function or not is also a question of interest. The 
probability of an abnormal behavior of this dependence at some Nb contents should 
not be excluded. If so, it would be useful to elucidate the reasons of such behavior. 
First-principles computer simulations of electronic structure and elastic properties 
of alloys is an important point making much more effective experimental works on 
a purposeful search for ternary alloys with needed elastic properties. The present 
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work targeted two goals: (i) first-principles calculations of electronic structure and 
Young's modulus of binary Tii_2,Nb2. /3-alloy at x=0.07-0.25, and (ii) determination 
of the partial substitution of Ti by Zr effect on electronic structure and elastic 
moduli of ternary Tii_x_yZTyNhx /3-alloy. 

Calculation detail 

The high-precision "ab-initio" FLAPW method realized in Wien2k package |3] 
was used for calculations of elastic constants Cu, C12 and Young's modulus of 
Tii_j.Nbj; /3-alloy. The exchange-correlation potential was calculated in generalized 
gradient approximation (GGA) according to the Perdew-Burke-Ernzerhof model 
[3]. Atomic sphere radii both for Ti and Nb atoms were 2.2 a.u.. The calculations 
were performed for 2,000 k-points in the first Brillouin zone. Inside the atomic 
spheres the wave function was expanded up to lmax=12. The electronic density 
and potential inside the spheres were expanded on the lattice harmonics basis up to 
Lmax=6. The RminxKmax parameter which controls the number of APW func- 
tions in the basic set equaled 8.35 that corresponded to 160 APW per atom. The 
Gmax value determining the number of plane waves in the inter-sphere potential 
expansion equaled 12. These parameters allowed to calculate the total energy of 
model ordered structures with accuracy of 0.001 eV. The calculations were car- 
ried out under conditions of a complete structural optimization considering both 
homogeneous and inhomogeneous deformations, i.e. the lattice parameters and 
atom locations were optimized. The model volumes were 3aox3aox3ao (54 atoms), 
where ao is the BCC lattice parameter of titanium. The elastic moduli were cal- 
culated for ordered structures Ti54_mNbm, (m=2,8,10,12,14) which corresponded to 
Tii_^Nb^(x=0.037, 0.148, 0.185, 0.222, 0.26) alloys. The large model volume used 
provides comparatively high accuracy in calculations of the atomic structure and 
deformations caused by substitution of Ti by Nb atoms in the alloys. The bulk 
compression, shear, and Young's moduli were calculated by the Voigt-Reuss-Hill 
method (VRH) [4]. For ab-initio simulation of elastic properties the approxima- 
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tion of an elastically isotropic material with C = where C = {Cu — Ci2)/2, 
was apphed. It should be noted that all calculations realized in the present work 
correspond to K. 

Results and discussion 

For a more adequate simulation of real alloys, Nb atoms in the model Ti54_mNbm 
structures were arranged as uniformly as was possible. Each Nb atom did not have 
any Nb atom adjoining. As mentioned above, the calculations of total energies 
for these structures were performed with a complete structural optimization. The 
components of elastic tensor for the cubic symmetry were found using the CubicE- 
last package (Wien2k |3]).The structures with tetragonal symmetry were generated 
from the model cubic structures for various Nb contents. Then appropriate de- 
formations with lattice parameters changes not higher than 3% were applied to 
these structures. The total energies were calculated by the FLAPW method for all 
generated structures. Using the total energies of initial and deformed structures, 
the components of elastic tensor for each structure were calculated. The calculated 
equilibrium lattice parameters and elastic constants of the /3-titanium and model 
Ti54_„Nb„ (n = 2, 4, 8, 10, 12, 14) are listed in Table [1] The concentration depen- 
dence of the Young's modulus, calculated for the binary /3-alloys, is shown in Fig. 
[H At low Nb contents negative values of the Young's modulus were obtained, that 
correlates with the phase field boundaries. At those low Nb concentrations initial 

Table 1: Lattice constants, elastic constants and Young modulus (GPa) of Ti (bcc) and Ti54_mNb„i(m = 2, 
4, 8, 10, 12, 14) structures modeling binary Ti-Nb alloys. 



Alloy 


Structure 


a, A 


Cn 


C12 


C 


B 


E 


Ti 


Ti 


3.2553 


79.80 


122.38 


-21.29 


108.19 


-68.35 


Tio.962Nbo.038 


Ti52Nb2 


3.2500 


104.58 


114.45 


-4.93 


111.16 


-15.02 


Tio.924Nbo.076 


Ti5oNb4 


3.2513 


138.31 


122.29 


16.02 


127.63 


37.68 


Tio.852Nbo.148 


Ti46Nb8 


3.2662 


145.71 


97.34 


24.19 


113.46 


67.74 


Tio.815Nbo.185 


Ti44Nbio(A) 


3.2586 


140.57 


97.16 


21.70 


111.63 


61.15 


Ti44Nbio(B) 


3.2715 


152.27 


94.10 


29.09 


113.49 


80.39 


Tio.778Nbo.222 


Ti42Nbi2 


3.2553 


146.67 


108.47 


19.10 


121.20 


54.43 


Tio.74lNbo.259 


Ti4oNbi4 


3.2553 


150.25 


108.84 


20.71 


122.64 


58.81 
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Figure 1: Experimental and theoretical concentrational dependences of Yung modulus of Tii-ajNb^;. 

stresses extensively deform the BCC lattice, allowing atoms to occupy energetically 
more favorable sites proper for the HCP phase. At higher Nb contents metastable 
martensite phases energetically more favorable than the a phase form in the real 
alloys. The calculations showed that the concentration dependence of the Young's 
modulus of the binary Tii_^Nb^ alloys is not a monotonic function. This depen- 
dence has certain peculiar features at Nb contents of x=0.15-0.18 (Fig. [T]). These 
results correlate well with the experimental data obtained in jU El El IB] and plotted 
in the same Figure. Up to date, the Young's modulus increase in the Nb concen- 
tration range mentioned above has been attributed to the presence of martensite 
phases. However, our results allow to suggest that the increase in Nb concentra- 
tion leads to considerable changes in the electronic structure of the ordered crystal 
structures, modeling alloy in question. It makes a significant influence upon the 
Young's modulus of the /3 phase. Another important factor affecting the Young's 
modulus is the distribution of Nb atoms in the model lattice. This can be demon- 
strated by calculating the Young's modulus for two types of Ti^^NbiQ structure A 
and B shown in Fig. [2l The A was used earlier for the calculations of the Young's 
modulus concentration dependence, see Fig. [H 

The B was constructed from the A by changing the Nb atom sites. In the B 
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Structure A 

NQ2T5 p.43m 



Structure B 

N^217 1-43 m 




#-Ti 
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Figure 2: Model Structures Ti44Nbio. 



each Nb atom has one or several Nb atom neighbors. The B structure appeared to 
be energetically more favorable and had higher Young's modulus as compared to 
the A, see the asterisk in Fig. [H It should be noted that the calculated concentra- 
tion dependence of the Young's modulus better correlates with experimental data 
if the modulus value for the B structure is considered (dashed line in Fig. [1]). The 
calculation results allow to formulate two important conclusions. First, the con- 
centration dependence of the Young's modulus of the binary single-phase /3-alloy is 
apparently not a monotonic function. So, the peculiar behavior of the experimental 
concentration dependence at certain Nb contents can be explained not only by a 
multi-phase state. Second, the distribution of Nb atoms in ordered structure mod- 
eling /3-titanium significantly affects the Young's modulus. The structures with a 
short range ordering, when each Nb atom has one or several Nb neighbors, have 
higher Young's modulus values as compared to the structures without such order- 
ing. To clarify the reasons of this effect, a detailed investigation on the nature of 
the chemical bonds in Ti-Nb system and the mechanism of the /3-phase stabilization 
by Nb is needed. Unfortunately, these problems are beyond the scope of present 
paper. The discrepancy between the calculated and experimental Young's mod- 
ulus values can be explained by using polycrystalline and multi-phase samples in 
the experimental works. Our results significantly differ from the theoretical results 
obtained in other works [9l [10], see Fig. [3] for comparison. Our theoretical values 
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Figure 3: Theoretical concentrational dependences of Young modulus of Tii-^^Nba; alloy in [9l [10] and this 
work. 

have the best correlation with experimental data, see Fig. [H The peak on the the- 
oretical concentration dependence at x=0.185 is a peculiar feature of our results. 
The better correlation with experiment is reached for the absolute values of the 
Young's modulus as well. The values obtained in j9l are significantly lower than 
both ours and experimental ones. The main reason of the better correlation of our 
results with experiment is possibly a larger model volume used (54 atoms), while 
in the works mentioned above model volume comprised 16 atoms (2aox2aox2ao, 
where ao is the BCC lattice parameter). In |j9j the FLAPW method for the Young's 
modulus calculations was applied as well. The elastic properties were simulated by 
the Voigt-Reuss-Hill approximation (VRH) [1] for polycrystals. In [10] the Young's 
modulus was calculated by the pseudo-potential method with the VASP package 
|11] in the elastically isotropic approximation. The Young's modulus of several 
model structures for fixed Nb contents but different atom distributions were also 
calculated in [TTj, the results are included in fig. [3l Unfortunately, the authors 
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did not describe the principles of structure selection used. The insufficient model 
volumes used in the two earlier works, in our opinion, is one of the reasons of a 
considerable discrepancy between the theoretical and experimental values. The 
three-fold larger model volume used in the present work apparently allowed a more 
adequate modeling the real atomic structure, and a more accurate calculation of 
deformation effects caused by substitution of Ti by Nb atoms. The concentration 
dependence of the Young's modulus of the binary Ti-Nb /3-alloy makes possible 
modeling the elastic properties of ternary alloys. Our aim is to predict the compo- 
sitions providing single-phase state and Young's modulus values lower as compared 
to the binary alloy. As mentioned in the Introduction, Zr as the third element is 
promising. Model ordered structures of ternary alloy were formed from the A and 
B structures shown in Fig. [2l The Nb content in these structures was constant, 
x=0.185. The Young's modulus calculated for this concentration was quite high. 
Table [H So, it is of interest to substitute some Ti by Zr atoms in proportions 
that would result both in a further reduction of the Young's modulus as compared 
to the binary alloy, and in single-phase state. Two model structures of the A type 
with differing Zr content, Ti44_mZrmNbio(m = 16, 28), were chosen. These struc- 
tures correspond to Tio.519Zro.296Nbo.i85 and Tio.296Zro.519Nbo.i85 alloys. As earlier 
for the binary alloys, the /3-type structures with a peculiar Nb distribution were 
constructed from the A structures. In the B structures Nb atoms had one or sev- 
eral Nb neighbors. As for the binary alloys, the calculations showed high energetic 
favorability of the /3-type structures. The calculations of the atomic structure and 
elastic constants were fulfilled for these four model structures. The results are 
listed in Table [2l The substitution of Ti by 30% Zr atoms in the Ti44Nbio(A) 

Table 2: Lattice constants, elastic constants and Young modulus (GPa) of Ti44_„iZr,„Nbio(m = 16, 28) 
structures modeling triple Ti-Zr-Nb alloys. 



Alloy 


Structure 


a, A 


Cii 


C12 


C 


B 


E 


Tio.519Zro.296Nbo.i85 


Ti28Zri6Nbio(A) 


3.3520 


130.00 


100.00 


15.00 


110.00 


43.04 


Ti28Zri6Nbio(B) 


3.3820 


144.22 


94.66 


24.78 


111.18 


69.19 


Tio.296Zro.519Nbo.i85 


Tii6Zr28Nbio(A) 


3.4253 


129.01 


94.67 


17.5 


105.67 


49.75 


Tii6Zr28Nbio(B) 


3.4315 


137.12 


89.70 


23.71 


105.5 


66.18 
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Figure 4: Region of existence of single phase and multi phase triple /3-phases Ti-Zr-Nb 

structure led to a reduction of the Young's modulus from 61.15 to 43.04 GPa in the 
Ti28Zri6Nbio(A) structure. Further increase in Zr content up to 52 at.% resulted in 
the Young's modulus of 49.75 GPa in the Tii6Zr28Nbio(A)) structure. Therefore, 
our results predict a significant reduction of the Young's modulus in the ternary 
alloys as compared to binary ones. However, the first-principles calculations can 
not predict whether the ternary alloys will be single-phase or not. To clarify this 
question we used experimental results of [12]. In this work the influence of Zr on 
the phase composition of ternary Ti-Zr-Nb alloys in a wide concentration range 
was investigated. The authors showed that partial substitution of Ti by Zr in bi- 
nary Ti-Nb alloy led to the /3-phase stabilization. The samples annealed at 1173 K 
were investigated. Based on the results of |12j, we plotted the dependence of the 
lowest Nb content providing single-phase state on Zr concentration, see Fig. |H 
This allowed to estimate the single-phase /3 field boundaries for the ternary alloys. 
The alloys modeled in the present work are designated by asterisk. The Young's 
modulus values calculated for the A-type structures are indicated. Apparently, 
the search for compositions providing low Young's modulus values is most effective 
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in the hatched area at moderate Nb contents. Experimental results of [ T3l [T4] are 
also included for comparison. The samples studied in those works comprised several 
phases, that resulted in comparatively high Young's modulus values. 

Conclusions 

The concentration dependence of the Young's modulus of several ordered 
Ti54_mNbm structures, modeling binary Tii_^Nba; alloys, was calculated by first- 
principles methods. It was shown that the features of this dependence at 0.15-0.18 
at.% Nb can be attributed to the types of Nb atoms distribution. A detailed 
comparison of the theoretical results with experimental data was fulfilled. The 
Young's modulus values were calculated for a set of ordered structures with dif- 
ferent Nb atoms distribution, modeling ternary alloys Tio.sis Zro.297Nbo.i85 and 
Tio.297Zro.5i8Nbo.i85- The results of these calculations allowed us to suggest the 
concentration region for single-phase ternary /3-phase alloys possessing low values 
of Young's modulus. 
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